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Abstract. We study the behaviour of the Standard map critical function 
in a neighbourhood of a fixed resonance, that is the scaling law at the fixed 
resonance. 

We prove that for the fundamental resonance the scaling law is linear. We 
show numerical evidence that for the other resonances p/q, q > 2, p 7^ and p 
and q relatively prime, the scaling law follows a power— law with exponent 1/q. 



1. Introduction. 

The standard map Jl5| , pl[ | is the area-preserving twist map of the cylinder: 
R / (2?rZ) x K given by: 

(ii) ( x '\ =t ( x \ = f x + y + esinx 

\y'J \yj V y + esinx 

where e is a real parameter. Note that for all values of e £ R one has SoT^oS = T_ e 
where S is the involution of the cylinder S(*) = ( X y)- This allows us to restrict 
the study of T c to e non negative. 

The dynamics of the Standard Map depends critically on the parameter e. For 
e = the curves yo = constant are preserved by the map and the motion on these 
curves is a rotation of frequency oj = yo- For e > not all these curves persist. 
Given an initial condition (xo,yo) we define the rotation number lot by: 

2-kujt (xo,y ,e) = lim 

71 — >oo n 

when the limit exists, where ni is the projection on the first coordinate: tt% (x, y) — 
x. For e > by Poincare-Birkhoff 's Theorem (see g ^ |T[ || and [Q for a 
different proof) homotopically non-trivial invariant curves with rational rotation 
number do not exist: the invariant curves of rotation number lot (x,y,0) = p/q, 
for integer p, q, q > relatively prime, split into 2q points when e > 0. This set 
of points is divided into two subsets of points each one composed of q points: the 
elliptic points and the hyperbolic ones. 
Let us define 

(■crit (w) = sup{e > : Ve' € [0, e] , there exists a C 1 , homotopically 

non-trivial, invariant circle of T e i whose rotation number is lu} 
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and e cr it (uj) = if u> G Q. The function w i— > e cr it (w) is called the critical function 
for the standard map. The definition of e cr u (uj) actually depends on the smoothness 
required for the invariant curve. However, whatever is the regularity of the curve 
(say C 1 , C k , C°° or C u ) this function is very irregular: it is almost everywhere 
discontinuous and zero on a dense G^-set. Moreover J.Mather proved p8| that 

su PweH e cnt (uj) < | and it is conjectured that sup wgR e crit (uj) = e crlt (^r 1 ) = 
0.971635... f2| |2|, §8 . 

For suitable irrational uj, KAM theory 24, ^, [l0| assures that for sufficiently 
small |e| there exists an analytic invariant circle upon which the dynamics of T e , 
restricted on this curve, is analytically conjugated to the rotation (9) = 9 + 2ttu>. 
This means that there exists an analytic (w.r.t. 9 on T = R/27rZ and |e| < eo, for 
some eo > 0) function u {9, to, e) such that 

x = 9 + u(9,uj,e), 1 + — (9,uj,e) > 



m 



V6»e 



and 9' = 9 + 2ttuj, namely u reparametrizes the invariant curve. Using perturbation 
theory we can expand this function into a power series w.r.t. e and Fourier' series 
w.r.t. 9 

n>l kGZ,\k\<n 

One may define the critical radius 



(1.2) 



which gives the so called Lindstedt series 
of convergence of such series as 



(1.3) 



Pcrit ( w ) 



sup lim sup|w„ (9, us) \ < 

9£T n^oo 



It is clear that p cr a (uj) < t C rit (uj), but it is an interesting open question to decide 
whether the two functions are equal or not on some set of frequencies. See j(| || for 
a positive numerical answer and Jl7| for a negative one. This will imply that the 
function defined for e € %< cr . it (w) as the sum of the series (L2) can be analytically 
continued outside its disk of convergence. 

In f3(| Marmi and Stark studied the relation between the t cr it (uj) and a "uni- 
versal function" depending only on the frequency: the Brjuno function B (uj) (see 



12 



|). Using the Greene residue criterion (||21|), Marmi and Stark find 



merically evidence that it exists (3 > and a positive constant C, such that 



(1.4) 



In e crit (uj) +(3 B (uj) 



< C 



uniformly in uj, as uj ranges in some set of quadratic irrationals, with B approxi- 
mately equal 0.9. Comparing their results with similar problems in |]35| , fi"7| they 
expected the value (3 = 2. 

In H and more recently in the behaviour of p C rit (uj) when uj goes to a rational 
value is studied; in both cases complex frequencie s we re used. It is well known that 
if the imaginary part of u is not zero the series fll.2| ) converges because there are 
no small divisors. In these papers the limit w — + 2 of p cr u (uj) when uj = 2 + irj for 
some real 77 — ^ was studied. The result is the following scaling law 



(1.5) 



Pcrit (uj) ~ Cp/ q 



for some positive constant C p / q depending on the resonance. 
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Finally very recently Berretti and Gentile (§] proved that 



(1.6) lnp cri t H+2BM 



for some positive constant C", uniformly w.r.t. ugR. This implies that the series 



given in (1.2) converges if and only if w is a Brjuno number (i.e. B (lo) < oo). 

The scaling law is related to this result which makes a link between rotation 
number and Brjuno function. One can prove ( |14[) that 



liminf [B(u;)+-\n 

cj is a Brjuno's number 



P 

LO 

q 



_ j 

- C p/q 



with ^ v j q some non negative finite constant depending on the involved resonance^], 
then the result of Berretti and Gentile Q implies the value (3 = 2. 

To end this short review about standard map critical functions and their scaling 



laws we cite the very recent work of Treschev and Zubelevich 46 , where they 
studied (quite general) area preserving twist maps of the standard cylinder into 
itself, e-close to an integrable ones. They proved that under some assumptions: 
given any rational number p/q £ Q, there exist two positive constants c and eo such 
that for all |e| < e , there exist two invariant homotopically non-trivial invariant 
C 1 curves with frequencies co t and w b and 

P 

toi < — < ujb and \loi, — lo[\ < ce. 

q 

The Theorem they prove (Theorem 3 p. 76 |Q) is based on a strong hypothesis 
(in particular assumption 2., which depends on the involved resonance) which is 
not verified by the standard map, except for the resonance 0/1. 

We consider the fundamental resonance 0/1, we fix e > and we define 

uj + (e) = inf{tj > : T e has a C , homotopically non-trivial, 

(1.7) invariant curve of rotation number to} 
and 

u>~ (e) = supjtj < : T e has a C , homotopically non-trivial, 

(1.8) invariant curve of rotation number to} 

note that lj + (e) = — lo~ (e). These definitions are similar to the one of given 
by R.S. Mac Kay (3^] in the study of the Hamiltonian of Escande. The existence of 
a stochastic layer around the hyperbolic fixed point x = y = 0, implies lo + (e) =/= 0, 
namely there is a "gap" in the curve which associates to each invariant circle its 
rotation number (frequency curve) . The main result of the first part of this paper 
is the following Theoremn. 



x Wc point out that the previous result is false if we replace liminf with lim. In fact considering 
the fundamental resonance for simplicity and taking the sequence of Brjuno's numbers (<^n) n >oi 
(jj^ 1 = a± (n) + n2 (^)_|_g i with a; : N — » N* strictly monotone increasing functions for i = 1,2 and 
£ any Brjuno's number. Then we have lim n __|_ 00 B (uj„) + lnaj n = if ln<J2 (n) = o (ai (n)) for 
n — > +oo, whereas if ai (n) = o fin ao (n)) for n — > +oo, then limn^+oo B (u> n ) + lnuj n = +oo. 

2 We point out that already in |34[ for a one degree of freedom Hamiltonian system depending 
periodically on time, logarithmic singularities were found for the breakup of invariant tori with 
frequencies closer and closer to rational ones. 
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Theorem 1.1. There exist two positive constant f2 /i and c /i such that for all 

\ne crlt (w) -ln|w| < c 0/1 . 



\u>\ < fio/i then 



Following p6| , to prove our Theorem we need a upper bound for In e cr a (uj) — ln| ui \ . 
We obtain this bound combining ideas taken from two different domains: the ex- 
ponential smallness of the splitting of the separatrices and the interpolation of 
diffcomorphisms by vector fields. 

The splitting of the separatrices for the standard map has been widely studied. 
Using the results of Lazutkin and Gelfreich jf{| |3f], |||, |33| we get a lower bound 
to the action corresponding to the last invariant torus. More precisely for a fixed e 
we define y + (e), respectively y~ (e), as the intersection point of the homotopically 
invariant C 1 curve of frequency ui + (e), respectively ui~ (e), with the j/-axis, note 
that y + (e) = — y~ (e). Then we prove (Proposition |2.2|) that 

2 

y + (e) > ce 2 ^ 

for some positive constant c and e small enough. 

Using a result due to Benettin and Giorgilli Q , we construct an Hamiltonian sys- 
tem whose time-1 flow interpolates the Standard map; this allows us to exhi bit an 
asymptotic development for the Standard map rotation number (Proposition 2. If ) 
which, joint with the bound on y + (e), leads to the desired upper bound. We note 
that to prove our result it will be sufficient to find an interpolating Hamiltonian 
system e^-close to the standard map; we nevertheless give the asymptotic develop- 
ment because to the best of our knowledge this result was not know and, secondly, 
the construction of an interpolating vector field exponentially close needs the same 
amount of difficulty as the construction of an e fc -close vector field, for any power k. 

In the second part (section_||) we present the Frequency Map Analysis (FMA) 
method of Laskar 1 18 , 2(], |25|, ^ and its application to the numerical investigation 
of the Standard map critical function. 

For the fundamental resonance 0/f , we find numerical results in agreement with 
the analytical ones presented in section |2[ For the other resonances we find the 
behaviour 



tcrit M — C. 



p/q 



P 

uj 

q 



for uj in a neighbourhood of p/q and c p / q > 0. We remark that using FMA we deal 
with real frequencies. 

Using the FMA tool we can do investigations directly on the critical function 
and go deeply in the resonant region. Actually we deal with frequencies that are 
only 6 • f CU 4 far from the resonant value (in |3(| the use of the Greene's residue 
criterion gives frequencies that are 10 -2 close to the resonant value). The result we 
show in this paper for the critical function is in agreement with the one of f36f . 

Supported by considerations on both methods (analytical and numerical) we 
make the following conjecture: 

Conjecture. For all rational p/q, we can find positive constant 
Cp/q, such that when uj is in a small real neighbourhood of p/q, then 
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Comparing this conjecture with ^ we see that the domain of analyticity with 
respect to e G C of the function defined by the series (1.2) is quite different from 



a circle as u> tends to a resonance It is an open question to understand 

completely its geometry: it could be an ellipse whose major semi-axis is in the 
direction of real e with a length proportional to the square root of the minor semi- 
axis, or probably a more complicated curve with the nearest singularity to the origin 
placed on the imaginary axis. We don't think our result and our conjecture are 
in disagreement with the one of J|, [| in fact we think there is a difference when 
dealing with real or complex frequencies. 

Acknowledgements. We are grateful to S. Marmi and D. Sauzin for the important 
help they provided during this work and to B. Chirikov and D. Shepelansky for 
some early discussions. Support from EC contract ERBCHRXCT94-0460 is also 
acknowledged. 



2. Scaling law for the fundamental resonance. 

The goal of this section is to prove our main result: Theorem |l.l| . We note 
that from the definitions of t cr it (u>) and u (e) we have e C rit (e)) = e then 
Theorem 



LI 



is equivalent to 



Theorem 2.1. Let uo^ (e) be the functions of e defined in (1.7) and (1.8). Then 
there exist e' > and c' > such that for all < e < e' we have 



lnlw* (e)| -Ine 



Using Theorem 3. page 76 of we can prove that there exist c > and eo > 
such that for all < e < eo there exist two invariant homotopically non-trivial C 1 
curves with frequencies lui and uib and 

bJi < < Wfc and \ujb — <jj{\ < ce. 



Recalling (1.7) and (L8) we get the lower bound 



lne - filial (e)| > In 2- lnc 



for e small enough. Thus the proof of Theorem 2.1 reduces to prove the following 
bound 



(2.1) 



lne-lnl^ (e)| < c" 



for some c" > and e small enough. This will be done in Proposition 2.4. Now we 
present the two results we need to prove Proposition [2.4| . 

Let us fix e > and recall the definitions of j/ ± (e): y + (e), respectively y~ (e), are 
the intersection point of the homotopically invariant C 1 curve of frequency u> + (e), 
respectively lo~ (e), with the y-axis. Then from the exponentially smallness of the 
splitting of the separatrices we deduce 

Proposition 2.2. There exist c > and e ' > such that for all < e < Cq we 

have 



(2.2) 



2v^r 
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We leave the proof of this Proposition to subsection 2.1 and we introduce the 
second element we need to prove our result. We rewrite the Standard map as 
perturbation of the identity map, namely 

, / sin x^ 



(2.3) 



z 
sin a; 



V 



we can pass from (2.3) to (1.1) putting // = yfe and [iz = y. And we define lo$ to 
be its rotation number, namely 



u>s (x, z, ji) = lim 



7T! (S»(X,Z)) 



when this limit exists. The relation between lo$ and lot is 
(2.4) co T (x,y,e) = y/ews (x,-^=,y/e 



Following M we construct (Proposition 2JS and Proposition 2.S ) an integrable 
Hamiltonian system whose time-l flow differs from the standard map in the formu- 
lation ( [2~3| ) for an exponentially small term: O (e~^^J for e — > + . This allows us 

to prove (Proposition 2.11 ) the existence of an asymptotic development of the stan- 
dard map rotation number in terms of y/e, whose first order term is the pendulum 
frequency for rotation orbits. In fact the best interpolating Hamiltonian system is 
an integrable perturbation of order O (y/e) of the pendulum Hq (x, z) = 4j- + cosx. 
This implies the following Lemma 

Lemma 2.3. Let ho be a real number strictly greater than 1. Let lot be the standard 
map rotation number and let Lo pen d be the frequency of a rotation orbit of energy ho 
for the pendulum H p (x, z) = 4- + cosx. For e > we consider the invariant circle 
of the standard map of initial conditions (xq, tjq) — (xq, -*^=^ , then for e sufficiently 
small we have: 



yo 



= o 



We can now prove the bound (|2.l|) 



Proposition 2.4 (Upper bound.). Let to ± (e) be the functions of e defined in (1.7) 
and (1.8). Then there exist e > and c" > such that for all < e < e we have 

lne-lnl^ (e)| < c". 



Proof. Lemma 2.3 and its refinement (Corollary 2.10) say that the standard map 
rotation number is "well approximated" by the pendulum frequency. Consider a 
rotation orbit of the pendulum starting at (0, zq). Its frequency is given by 

Wpend (^o) = 77" 



where k 2 = 



(i + 4)- 



hn = 1 + % and 



kfC{k) 



dip 



yj\ — k 2 sin 2 i\) 
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is the complete Elliptic integral of first kind. For k = 1 (that is zq = 0) the integral 
diverges logarithmically and (see fij], formula 17.3.26 p. 591) 



lim 

zo^O 



JC-ilni 



= 2 In 2. 



Therefore in the limit zq 

(2.5) Upend 



0, Wpend behaves as 



In- 



= 



ln|z 



Assu me n ow that ur (0, yo, e) exists for e small enough, where yo — zo/i. By (2.4), 
Lemma 2.3 and (2.5) one has 



(2.6) 



I^t (0,y o ,e)| > 



In 



8V7 



O 



(H«) 



If yo is the intersection of the last invariant torus wi th t he axis x 
yo = y + (e) and lot (0, yo, e ) = w + (e), then the bound (2.2) gives 

n 



0, namely 



\u>+ (e)\>Ve 



C 



line 







¥) 



for some positive constant C. Then for e small enough we conclude that there exist 
another positive constant c such that 



k + (e)|>ce. 

Recalling that lo~ (e) = — lu + (e), the Proposition is proved. 



□ 



2.1. Splitting of the separatrices. The aim of this subsection is to introduce 
some results about the splitting of separatrices for the Standard Map necessary to 
prove Proposition 2.2. The splitting of the separatrices has been studied by many 
authors, here we use the notations and the presentation given in |l9| |. 

For e = the circle y — is formed by fixed points of the Standard Map, for 
e > only two fixed points survive: (0,0) and (tt, 0). The first one is hyperbolic 
(for these values of e) and the other one elliptic. The linearised map at the origin 
is 



1 + e 

e 



and its eigenvalues are A and A 1 where 

A= 1+^+ 



The stable, W s , and unstable, W u , manifolds of this fixed point are analytic 
curves passing through (0,0). We introduce the parameter 6 = In A and we note 
that e ~ S 2 for small e. We parametrise the part of the unstable, W" separatrix 
which grows upwards from the origin, with (x~ (t) , y~ (£)), imposing the boundary 
conditions^ 

lim x~ (t) = 0, x" (0) = 7T. 

t— * — GO 



3 The y (t) can be found using y (t + 5) = x (t + S) — x (t). 
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We assume that t — is the hrst intersection of W" with the vertical line x = n. 
The parametrisation of Wf is given by 

(x + (t) , y + (t)) = (2tt - a; - (-t) , y" (-t) + e sin a;" (-*)) . 

One can show that 



lim x 

t — >+oo 



(t) = 0, a;+(0)=7r, 



i.e. t = corresponds to a homoclinic point. Lazutkin 32 proposed to study the 
homoclinic invariant defined by 



fl = det 



±- (0) x+ (0) 
y-(0) y+(0) 



The homoclinic invariant has the same value for all points of one homoclinic trajec- 
tory and it is invariant with respect to symplectic coordinate changes. Gelfreich in 



19 1 , Proposition 2.1, gives an asymptotic expansion in power of 6 of the parametri- 
sation x (t). Its first terms are 



X (t) + Xl (t) + X 2 (t) 



(2.7) =4arctane* 



, 1 sinh t 
4 (cosh tf 



41 sinht 



91 sinht 



1728 (cosh tf 864 ( cos h t) 4 



Of course the zero order term is nothing but the parametrisation of the separatrix 
of the pendulum. 




Figure 1. The stable and unstable manifolds for the Standard 
Map with some of their infinitely many intersection points. 



The main result of |19| is f2 > jie~~ for some positive constant c and 8 small 
enough. Using this result we obtain an estimation of y (e) in Proposition 2.2. 
Proposition 2.2. There exist c > and Cq > such that for all e < 6q we have 



ll/* («)l > ce~ 



Remark 2.5. Compare with Figure |71| where we plot ln|y + (e cr it)—y 
£crit, with data obtained numerically using the FMA. 



{e cr u)\ against 
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Proof. For a fixed value of e > 0, the stable and unstable manifolds intersect firstly 
at (tt, y (e)) and then infinitely many times. We call lobe the arc of manifold between 
two successive intersections and lobe area the area bounded by the lobe. Each lobe 
will be distinguished by an integer: the number of iterates (of T e ) needed to put 
the point (n, y) onto the right point base of the lobe. 

The hyperbolic fixed point (0, 0) stretches the unstable manifold along the di- 
lating direction (corresponding to the eigenvalue A), so there exists an integer k 
such that the k-th lobe intersect the y-eods at some point y, n t (e). Because the 
invariant circles of the Standard Map are at least Lipschitz graph, we cannot have 
homotopically non-trivial invariant circles, passing by (0, y') with y' < yi nt (e), thus 
y + (e) > yi n t (e). The lobe area is conserved under iteration of the Standard Map 
and by |l9), Corollary 1.3, this area isQfor 6 small enough 

Alobe ~ e T 

where lj = 1118.827706 • • • and S = In A ~ y/e for e small. 

In Figure || we represent the geometry of the generic fc-th lobe: the points 
Ak = (x-k,y~k) and Bk = (-^-(fe-i) , 2/-(fe-i)) are the iterates of the standard map 

sc-(k-i) = x-k + y-(k-i) 

f-(fc-i) = V-k + <5 2 sina;_ fe 

the segment of line A^Ck is tangent to the unstable manifold at Ak, while AkBk is 
tangent to the stable manifold at Ak- 




FIGURE 2. Two intersection points of the stable, W s , and unsta- 
ble, W u (dotted curve), manifolds, producing the fc-th lobe. 



4 In the follow we will use often the standard notation (see PQ ] n.4 chap. II) / (i) ~ g (t) for 
t — » 0, this means that / and g are equivalent when t is close enough to 0. We can restate it by 
saying: f (t) = g (t) + o(g (i)) when t -> 0. 
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The angle ctk is defined by 

y-(k-i) -V-k _ <5 2 sinx_ fc 



tan at; 



-O-i) - x-k 



An 



Let us consider the lobe which intersects for the first time the y-axis and let 
us call yi n t the intersection point. We don't know exactly the analytical form of 
the arc of unstable manifold joining At = \x_z,y_£) to (0, y% n t)i so we can not 
calculate exactly yint- Then we approximate this arc with its tangent at A k : 

(2.8) y = y_ k +(tanp k )(x-x_ k ). 

This straight line intersect the y-axis in some point y and |y — yi n t \ = O I x 2 



kj 

We approximate the lobe with a parallelogram of sides AkBk and A^Ck where: 

1. the straight line through Ak and Ck is tangent to the lobe at Ak] 

2. the area, AkBk ■ A k C k sin ((3k — £*fe), of the parallelogram is equal to Ai b e . 



The condition of intersection with the y-axis is now A k K k = x_ k : this defines 
implicitly k. 

Let Hk denote the orthogonal projection of Ck on the straight line through Ak 
and B k . Then A lohe = A k B k ■ C k H k = ^^C k A k sin {(3 k - a k ). Since A k K k = 
— CkAk cos/3fc, the intersection condition becomes 

— — cos at cos 8t 

2.9 A k K k = - . ) a k s Aobe = x_ k . 



From (p.Tj), using t = —kS one gets 



thus 

tano-g ~ 6 (l + e~ 2kS ^J , sin a k ~ 5 (l + e~ 2kS ^j , cosa fe ~ 1. 

To estimate [3 k we use again the fact that the intersection will happen very close to 
the origin: the tangent to the unstable manifold at A k is very close to the expanding 
direction at the origin, so we get 



tan/% ~ -5 ( 1 + £ ) , cos/% ~1-<S 2 . 



Using these approximations we can rewrite (2.9) as 

e 

from which we obtain 

e2 -2fc(5 



(1 _ s 2 ) ^e-^ = Ae~ U ASe k \ 26 + o (5 2 e~ 2kS ) 

V ' 7T l+e -2k8 \ J 



for some positive constant C±. This allows to compute k explicitly. To determine 



y, we recall (2.8) and using the previous value of k we get 



.'/ = //_, ■'• , ian /, - iSe~ ks + Ae-' kS S ( 1 + ^ )- C\< 
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for some new positive constant C[. Finally, since \y— yi n t\ = O (x 2 
we conclude 



= 0[ e~7? 



-V (e) =y + (e) > ce ^5, 
for some positive constant c. □ 

2.2. Interpolating vector fields. Once we have determined a lower bound for 
y + (e) in terms of the intersection of the lobe with the y-axis, we need to compute 
a lower bound for the rotation number of the corresponding invariant circle. To 
this purpose following Q, we consider the problem of the construction of a vector 
field interpolating a given diffeomorphism. With this subsection and the following 
we prove the existence of an interpolating Hamiltonian system exponentially close 



to the ma p (|2.3| ) (Proposition 2J3 and Proposition 2.9). These results will be used 



in section 2.4 to prove the existence of an asymptotic development for the stan- 
dard map rotation number (Proposition 2.11) from which Lemma 2.3 immediately 
follows. 

Let I? be a bounded subset of R" 1 for m £ N* , la. real interval containing the 
origin and let us consider the one-parameter family of diffeomorphisms $ 

(x,fi) i-» f {x,n) 
such that for all x € D and fi € I it has the form 
(2.10) ^(x, fJ ,) = x + J2^n(x)fJ, n . 

n>l 

The problem consists in finding a one-parameter family of F : D x I — > IR m 
vector fields such that the time-1 flow of the differential equation 

I -"<*■"> 

coincides with the diffeomorphism 



We assume that ^ n (x), defined in (2.10), are real analytic functions of a: in some 



complex neighbourhood V p of D, defined as follows. Given p = (pi,... ,p m ) S R™ 
then 

V p = (J {z e C m : \zi - x z \ < Pt ,l<i< m}. 
xeD 

We introduce standard norms on complex analytic functions. Let / be a scalar 
analytic function on some domain D c R m , let p £ R™ and consider its complex 
extension T> pi then 

(2-11) = sup |/(z)| 

z£T> p 

defines a norm. If F is vector valued analytic function F = (fx, . . . , f m ) on T> p , for 
some p = (pi, . . . , p m ) € R+, then we define a norm (noted as in the scalar case) 

||F|| P = max \\fiWp. 

l<«<m 

If is a formal /Lt-series we can look for a formal interpolating vector field 
F(x,p) = ^2 n>1 F n (x) p n whose time-1 flow coincides with the diffeomorphism, 
as we show with the following Proposition. 
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Proposition 2.6. Let (x./i) be a diffeomorphism of the form (2.10) and assume 
^ n (x) be analytic in T> p for all n > 1 and for some p G i?"' . Then there exists a 
formal vector field F^ of the form 



(2.12) 

such that 
(2.13) 

for all x £ T) p . 



e L *v ( x ) = * (x, fj,) 



Proof. The proof is a very natural application of the Lie series. We give it for the 
sake of completeness. 

By definition of exponential of the Lie derivative 



n>l 



Looking for in the form ( [2.12| ) and using the linearity of the Lie operator we 
have 



n>l 

For all integer n > 1 and 2 < m < n we define 



n— m+1 



£i,n (x) = L Fn (x) = F n (x) £ m , n (x) = ^ L Fj L m -i tn -j{x) 

3=1 

and it is easy to prove that 

n — 2-^i " •'-m.n- 

n>m 

Then by (|]l|) and ( fc.lO| ) we conclude that for all k > 1 

V -:£„,fe (at) = (x) , 

namely 

(2.14) L Fk (x) = t> k (x) - V — £ n , fe (x) . 



Let 



*r'J and let G = £* =2 ^n.fc (x) = (d, . . . , G m ). Note 
that G depends on ^ j and i^-i for j < fc. Using ( 2.14 ) we have for all 1 < Z < m 

FlP {x)=*f+G h 



Am) 



n=2 



,fc i 



from which we obtain 



inF fe (x) = (F« (x),...,Ff } (x)). 



□ 



If ^> is analytic w.r.t /i then we have the following result ( B Proposition 1 p. 1122) 
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Theorem 2.7. Let ^(x,fi) be a diffeomorphism of the form (2.10) with \& n (x) 
real analytic in some complex domain T> p for p G K."', and let F^ be the formal 
interpolating vector field obtained in Proposition 2.t. Assume that for all k 

ll* fc ||p<r 7 fc -\ 

for some positive constants T and 7 then 

N^xlU < r, vfc>2 ||Ffc||| < ±k k - 1 k - 1 r, 

with (3 = 4max{7,r}. 

Moreover if \& (x, ji) is symplectic then for all n F n (x) are locally Hamiltonian. 



In the form (2.3) the Standard map verifies the hypothesis of Theorem 2.7, but 
to get optimal estimates we state an Hamiltonian version of Theorem 2.7 adapted 
to the Standard map. 



Proposition 2.8. Let z > and consider the interval L 
and Sq > we define the complex cylinder 



-- (—z, z). Fixed Rq > 

V ( Ro,S ) ={x£ C/2ttZ : \$Sx\ < S } x (J {z 1 e C : \z - z'\ < R, }. 

zei 

Then we can find a formal Hamiltonian H (x,z,fj,) = YlnLi Hn-i O^- 2 )/- 4 " such 
that 

V (x, z) e V (R ,s ) e LH ( x , z ) = Sn ( x , z ) ■ 

Moreover for all < d < 1 there exist A = 2So\\H \\ {RatSo} and B = f " H °jj^' So) 
such that 

Vfc \\H k \\ {1 „ dKRo , So) <AB k (k + l)\, 
namely H is 1-Gevrey w.r.t. fi. 



With norm defined in |2.11| it results that: ||i?o||(_R ,s ) = + e °. 
The proof of Proposition 2i is an adaptation of the proof of Theorem [2.7| and 
we omit it. 

We write down the first few terms of H 

(2-15) 



H (x,z, n) = — + cos x + ix— sin x + O (/1 2 ) 



and as we stated this is a perturbation of order O (y/e) of the pendulum. 

2.3. Exponentially small estimates. Proposition gives us an Hamiltonian 
function represented as a (maybe divergent) series. It is then interesting to define 
a truncation of this series, namely take N S N* define ffW (x, z, fx) by 

N 

ffW (a:, z,/i) = ^ H n _i 

n=l 

and to ask for the relations between 5 M (x, z) and vj/t^) (x, z, fi), where the latter is 
defined by 

W (N) (x,z,fi) = e L H(») {x,z) 

We can prove that we can chose N such that S l _ l (x,z) and vpW (x,z,(i) are 
exponentially close. More exactly: 
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Proposition 2.9. Let i?W (x, z, fx) be the truncation of the 1 -Gevrey interpolat- 
ing Hamiltonian system given by Proposition |<?.<ji| , and let (x, z, fi) be its time-1 
flow. Then we can find a positive integer N* such that for all < d < | and for 

all | M | < 2^ 

||5 M (a;^)-^ JV *)(a ; ,^M)||( 1 _| d ) ( _ R0iSD) < |/i|i*rw, 

with D = M + e ^j;™ d i , D' = (2eS) _1 , M = max{i? + |^|e 5 °,e So } and to = 
max{i? , e s "}. 

This Proposition is an adaptation of the corollary 1 p. 1124 of [Q and we don't 
prove it. 

We will note in the follow H* (x,z,/i) = H^ N > (x,z,fi): the "best" truncation 
of the interpolating Hamiltonian. 

Let us introduce the frequency of the rotation orbit of H* . 

TTi (e" £ «- {x,z))-x 
2iru>* [x, z, /J,) = lim . 



Then from Proposition 2.9 we deduce the following corollary. 

Corollary 2.10. Let u>s be the rotation number of the map (2.3). Let H* be best 
truncation of the interpolating Hamiltonian system given by Proposition and let 
27rw* be its frequency. Then we have 

(2.16) \w s (x,z,n)-u m (x,z,n)\ <D\n\e-^. 

Proof. We can prove that for all positive k and for initial data (x, z) on some domain 
such that we can iterate e L "* and at least k times, we have 

fe-i 

(x, z) - e kL »' (x, z)=Y, e (k - l - 1)L «' (S, - e L »' ) S< (x, z) . 

1=0 

For real x, z and /x, the standard map is real valued so = S„ (x, z) is in 27rT x R 
for all I > 0. 



Let us define % = (S M - c Lh ') (fj), then Proposition pT9| implies < D\^\e ra 
for (x, z) € T^(i- a d)(Ro S )- Such a D depends on the domain troughs d, for some 
fixed So and Rq, but using the previous remark on the reality of £/, we can fix once 
for all a domain, and so a constant D, in such a way that £j belongs for all Z to this 

domain. For all n: |e TlL «* < D\[i\e~W\ 1 because e" L «* ($) = + O (#). 
We conclude that 

[tti (5* (a;, z)) - z] - [tti (e fei ** (x, z)) - i] | <k\n\De-&, 

from which @.16| ) follows. □ 

2.4. Asymptotic development. We can now prove the existence of an asymp- 
totic development for the standard map rotation number in a neighbourhood of 
e = 0. 

The Hamiltonian H* defines an integrable (it has only one degree of freedom) 
Hamiltonian system whose period is given by 

27r f 27r dx 
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where ho is an energy level corresponding to rotation orbits. 

Let us fix ho > 1. We are interested in finding curves on the energy level 
H* (x,z,fi) = ho- Take zo > and consider the point Pq = (xo,Zo,0) such that 
Po £ H*~ x {ho}, from Q2.15D it follows that: ^ (P ) = z > 0. We can then apply 
the implicit function Theorem which assures that it exists a function z = z (x, fx) 
such that locally it satisfies 

(2.18) H* (x, z, n) = h Q 

moreover this function is smooth. We compute the first terms and we get: z (x, fi) = 
y / 27^o^-^os~r) — ^ sin a; + O (\^\ 2 ) and 

dH* 



dz 
(2.19) 



(x, z, fi) — y2 (ho — cos a;) 



2 



N* 



Q IT 

sinx + J2 ^JT (x ' f) " ^2(^o- cos a;) = O (|/x| 2 ) 



n=2 







Recalling 2.4 we prove the following Proposition 

Proposition 2.11. Let LUx(x,y,e) be the rotation number of the standard map 
in the formulation (1.1), then it exists a formal y/e-power series X^^Lo (V*) 
asymptotic to ujt (x, y, e) for e — > 0. Namely for all positive integer N we have 



(2.20) 



lim (Vi) 



-(iV+1) 



A? 



E B « (^) 



n+1 



Proof. The solution, z, of (2.18) given by the implicit function Theorem is analytic 
w.r.t x € T and \i in a neighbourhood of 0, AP- (x,0) is positive for ho > 1, then 



for /it close enough to the origin 
and /i. We can then write the following Taylor development 

-l 



(x,fj,) is positive and analytic w.r.t to x G T 



(2.21) 



dz 



(x, fj) 



^K n (x,h )^ n , 



n>0 



where (K n ) n>0 are some known functions of (x,ho). Using (2.19) we get 



(2.22) 



K (x, h Q ) 



y2 (ho — cos x) 



and ifi (x, /iq) = 0, 



we left to remark 2.12 an iterative schema to calculate K„ 



For all k > we can integrate the K n (x,ho) and we rewrite (2.17) as 

— = 2 7 rVC7„(/ l0 )M", 



from (2.22) we get 
Co (h ) 



1 



and d (ho) = 



OJpend (ho) Jo \/2 (ho — cos x) 

where uj pen d (ho) is the frequency of a rotation orbit of energy ho of the pendulum 
H p (x, z) = 4j- + cosx. 
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Clearly for all positive integer N we have 

JV 



(2.23) 



lim u 



-JV 



1 



w* (h 



n=0 



= 0. 



We show now that we can obtain a similar estimate for cj*. Let us set, for all 



positive integer N, Sn (a*, ho) = SnLo ^« (^o) A*™- Because Sjy (0, /io) 7^ we can 
express 1/Sn as a [i-power series 



1 



1 



E^" E 



(2.24) =5^B B (fto)/i", 



^ C ni (h ) . . . C nm (h ) 

-+n m = 
2<m<N 



n>0 



where the second line is the definition of the B n 's, Bo (ho) — Co(ho) 1 and 
B 1 (ho) = 0. Recalling ( p723| ) and ( gTgg ) we can then conclude that for all pos- 
itive integer TV 



JY 



/ n>JV+l 



n=0 



SW (a*, fto) 



o(\ii 



and 



n=0 



= 0. 



From |2.4| and corollary [2.10] we get for all positive integer N 

N 



^ - E B » ^0) W 



n+1 



n=0 



< eDe 



JV 



dividing by (\fe) N+1 and passing to the limit f^flwc obtain (2.20). 



□ 



From this Proposition Lemma 2.3 follows immediately. 



Remark 2.12. The calculation of the K n 's defined by (2.21) is a little cumbersome 
but easy, it needs only some algebraic manipulations of power series. We assume 
known the solution, z(x,fj,,ho), of ( 2.1 8| ) given by the implicit function Theorem, 
and using the analyticity hypothesis we write it as a fi-power series 

z (x, fi, ho) = zo (x, h ) + 2J Zn (x, ho) 

n>l 



Starting from H* given by (2.15), we calculate its partial derivative w.r.t. z, 
then we substitute z = z and we reorder the powers of /i to get 

(x,z,fi) ^ zo(x,h Q ) + 2_^ z ilJ- + Z^l 1 Z_^Ti — — — ( x > z o) 



dz 



l>2 



E 



l>3 



Z0 



n-\-m—l k—1 
i>l,2<n<N 



k\ dz 



(x, h ) + ^Ti-i (x, h ) /J, 1 . 



mi + ...+m k =m 
m,>l 
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This series is invertible in the field of the formal fi-power series, then (2.21) defines 
correctly a \i-power series, which is given by 

«£ (M ,„ ) )- 1 _i +E ,,£fc£ E H Hu 

' u l>2 m=l ii+...+/ m =i 



2i>2 



= if (i, /i ) + ifj (x, /i ) M ! 



i>2 

3. Numerical Analysis of the scaling law at the resonance p/q. 

3.1. Frequency Map Analysis. The Frequency Map Analysis of Laskar (FMA) 
is a numerical method which allows to obtain a global view of the dynamics of 
Hamiltonian systems by studying the properties of the frequency map, numerically 
defined from the action like variables to the frequency space using adapted Fourier 
techniques. Thanks to its precision it was used for the study of stability questions 
and/or diffusion properties of a large class of dynamical systems : Solar System 
|26|, Particles Accelerator fi(J, Galactic dynamics @ El, Standard Map j|| |§. 
We present here the outlines of the method, following jg7). 

3.1.1. Frequency Maps. Considering a n-degrees of freedom quasi-integrable Hamil- 
tonian system in the form 

(3.1) H (J, 9; e) = H (J) + eHi (J, 9) 

where H is real analytic for (J, 6) G B x T", being B an open domain in K", and 
e is a real "small" parameter. For e = this system reduces to an integrable one. 
The motion takes place on invariant tori Jj = Jj (0) described at constant velocity 

v i (J) — in 2 - , for j = 1, .... n. Assuming a non-degenerate condition on Hq 

oj j J(0) 

the Frequency Map F : B -> R" 

F : F( J) = v 

is a diffeomorphism on its image f2. In this case KAM theory |24|, |[ [hi) insures 
that for sufficiently small values of e, there exists a Cantor set f2 e C f2 of frequency 
vectors satisfying a diophantine condition 

C 

l< k,v>\ > 



for some positive constants C e and m, for which the quasi-integrable system ( |3.lD 
still possess smooth invariant tori, e-close to the tori of the umpertubed system, 
with linear flow t \—> Vjt + 9j (0) mod 27r for j = 1, . . . , n. Moreover, according to 
Poschel Q there exists a diffeomorphism * : T™ x ft — > T" x B 

*:(^, I /)^(^J) 

which is analytic with respect to <f> in T™ and C°° w.r.t v in £! e , and which transforms 



the Hamiltonian system (3.1) into 
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For frequency vectors z/ G f2 e the invariant torus can be represented in the complex 
variables (zj — Jje lSj )j = i^ n by a quasiperiodic expression 

(3.2) Zj (t) = Zj (0) e^ 4 + Hm M e l< "^ >4 

m 

If we take the section 6 = 6*n, for some 8q E T™, of the phase space, we obtain a 
frequency map Fe : -B — > £1 

(3.3) F e „ : J^^^- 1 ^, J)) 

where 7T2 (i/>, f) = v is the projection on Q. For sufficiently small e the non- 
degeneracy condition ensure that Fg is a smooth diffeomorphism. 

3.1.2. Quasiperiodic approximations. If we have the numerical values of a complex 
signal over a finite time span [— T, T] and we think that it has a quasiperiodic 
structure represented by a quasiperiodic function / (t) , we can construct a quasi 
periodic approximation, /' (i), as follows. We represent the given signal as: 

(3.4) f( t ) = e i ^ t + J2 a fe e l< ^>* a k e C, 

feez«\(i,o,...,o) 

and the signal reconstructed by FMA as /' (<) = ^fc=i CL k e tuJkt . Frequencies u>' k and 
amplitudes a' fc are determined with an iterative scheme. We first determine the first 
frequency to^ as the value of a which maximise: 



4>{<^) = \^ f_ T f{t)e- M X {t)dt 



where \ (t) is some weight function^]. Once the first frequency is found the associ- 
ated complex amplitude a\ is obtained by orthogonal projection. We then iterate 
the scheme to the fu ncti on: fx (t) = f (t) — a' l e lu ' lt . 

Given the signal ( |3.4| ) over a time span [— T, T], if we assume some arithmetic 
properties of the frequencies (Diophantine numbers for example), then one can 
prove (p^|) that, using the Hanning filter of order p: \p (t) — 2 (2p)! + cos7rf) p , 
FMA converges towards the first frequency with the asymptotic expression for 
T — > +oo: 



A T 2 p+ 2 ^ n 2p+1 \T 2 P-t 

k k 

where Clk —< k, v > —v\ and A p = [\ — J2k=i w\ We can t ^ lus notc 
that the accuracy of the determination of the frequency is 0(l/T 2p+2 ), which will 
usually be several order of magnitude better that with usual FFT ( 0(\/T)). 



3 A weight function is a positive and even function such that: 

cT 



In all the computations we will take the Hanning filter: \i (*) = ! + cos 9- 
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3.1.3. Frequency Map Analysis. With the previous algorithm, it is possible to con- 
struct numerically a frequency map in the following way : a) Fix all values of the 
angles 9j = 9jo- b) for any value Jo of the action variables, integrate numerically 
the trajectories of initial condition (Jo, Go) over the time span T. c) Search for a 
quasiperiodic approximation of the trajectory with the algorithm of the previous 
section, d) Identify the fundamental frequencies v of this quasiperiodic approxima- 
tion. The frequency map is then 



(3.6) Fl-.J^v 

and from the previous section, we know that for T — > +oo, on the set of regular 
KAM solutions, Fj Q — > Fg Q . In particular, Fg Q should converge towards a smooth 
function of the set of invariant KAM curves. The destruction of these invariant 
curves can thus be identified by the non regularity of the frequency map Fj Q . 

3.2. Application to the Standard Map. We are interested in the study of 
homotopically non-trivial invariant orbits of the Standard map, in particular in the 
relation between the rotation number w and the critical function e cr n (lu). Given 
a couple of initial conditions (xo,yo) the FMA gives us approximate amplitudes 
and frequencies of the resulting orbit. For e = all invariant curves are tranverse 
to section x = constant, for e sufficiently small this fact still holds, that is the 
tranversality property. For e > 0, "curves" of rational frequency do not persit and 
are replaced by "set of periodic fixed points" such that if u(xo,yo) ~ p/q then 
this set contains 2q periodic points, an half of which are hyperbolic (HFP) and the 
other elliptic (EFP). 

Near a HFP appears a stochastic layer, that is there are not invariant curves 
with frequency "very close" to the frequency of the HFP. This is shown in Figure |^ 
where the stochactic layer is revealed by the nonregularity of the frequency map 
around the frequency value 1/3. 



For a fixed value of e and a HFP of rotation number p/q, we study the frequency 
curve y i— > lu (x, y) for y in some neighbourhood, X, of yapp, and x = xhfp fixed. 
We generalise in an obvious way the definitions (1.7) and (1.8) of Lu + (e) and cu~(e) 



defined in the vicinity of rotation number for invariant curves in the vicinity of a 
general HPF of rotation number p/q as w+ /? (e), w~ /g (e) and Auj p / q = uj+ /q - 

Similarly we extend the definitions of y ± (e) as y p , q (e). 

Numerically y p / q (e) is determined as the smallest value of y(0), larger than yHFP, 
for which the frequency curve is still regular (see appendix [A]) , and similarly for 
Vp/q ( £ )- Tlicn we assume uj± /q (e) = lu {x HFP ,y± /q (e)) . Studying Auj p/q (e) as a 

function of e we are then able to reconstruct the curve e cr a (Aoj p / q ) for a fixed 
resonance. In fact the non-existence of invariant curves, for a fixed value of e, with 
frequency lu, between and tu^ q means that e cr u (lu) > e and e cr u i^^/^j = e - 

3.3. Presentation of the numerical results. Using the FMA we studied the 
behaviour of e cr ;t (ALu p / q ) in neighbourhood of the resonances f>|>|>§>3>3- The 
results are presented in Table |l| and in Figure ||, where the plot of In e cr it (Aw p /,) 
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Figure 3. Particular of a nearby zone of a HFP, section x = 4.408 
and e = 0.7 (neighbourhood of the 1/3 resonance). Note that for 
y G [2.535, 2.545] the frequency curve looks very irregular. 

as a function of In Aw p / 9 is shown. All numerical calculations are done with 100 
digits accuracy using the Fortran 90 multi-precision package mpfun ||. 



These results suggest the following geometry of the t cr i t -io plane. Consider only 
the resonance p/q, then for to close to p/q we have: 



(3.7) 



tcrit (w) 



C 'p/g 



P 

u> 

q 



Because of the density of the rational numbers on the real line, close to p/q there 
are infinitely many rational numbers p'/q', with q' bigger and bigger as p'/q' is 
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i.O -6.0 



-2.0 



Figure 4. Plot of lne crit vs In (Aw p/(? ) for p/q e {0/1,1/5,1/4,2/5,1/3,1/2}. 



closer to p/q. The relation between e cr it and u> is not the simple one given by (3.7), 
even in a very small neighbourhood of p/q, but we must add the contributions of 
all p' /q' near p/q. These contributions are nevertheless negligible because for very 
big q' the right hand side of (3.7) is almost equal to one. Then the local description 
given by the (3.7) should be a very good approximation of the relation between 
e cr it and to (see Figure pi). 



4. Conclusions. 



We have studied both analytically and numerically the Standard Map critical 
function in a neighbourhood of a fixed resonance. For the fundamental resonance 
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resonance p/q 




Cp/q 


0/1 


0.937 


1.066 


1/5 


0.175 


0.489 


1/4 


0.223 


0.723 


2/5 


0.175 


0.714 


1/3 


0.312 


1.000 


1/2 


0.427 


0.957 



Table 1. lne cr ;t = A p / q hi/S.uj p / q + C p / q , where the coeffcients 
A p / q and C p / q are determined by a linear least square fit to the 
FMA results. 




we find a linear scaling law, whereas in Q and || a quadratic scaling law was 
found. We think these two results agree, in fact we think there is a difference in 
the scaling law for real or complex frequencies. This shows in particular that for 
diophantine values of the rotation number to, sufficiently close to zero, the critical 
function e cr it{u>) and the critical radius of convergence p C rit{^>) are different. More 
precisely e crit {uj) > p crii (w) > 

We stress again that this result implies that the conj ugating function u defined 



as the sum of the Standard Map Lindstedt series ( |l.2|) for |e| < p cr it{^>), can be 
analytically continued outside its disk of convergence at least for e real. This implies 
that in the plane e-complex the domain of convergence of the Lindstedt series is 
not a disk; it will be interesting to study its geometry, in particular its boundary: 
it could be a fractal curve as in the case of the Schroder-Siegel center problem jl3| , 
we only know that the nearest singularity to the origin is on the imaginary axis. 



The proof of Theorem 1.1 is different from the one of §. The FMA method 



doesn't allow us (at least for the moment) to use complex frequencies so we couldn't 
verify the results of directly. It would be interesting to investigate the scaling 
law for the generic resonance p/q with our methods, but the proof of Theorem [0] 
is specific to the fundamental resonance: we don't have a bound of the initial data 
for the last invariant torus for a generic resonance as we proved in Proposition pT3 
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for the resonance 0/1. Supported by the numerical results of section ^, we can 
nevertheless make the following conjecture 

Conjecture 1. For all rational p/q, we can find a positive constant c p / q , such that 
for all lo in a small real neighbourhood of p/q 



(4.1) e cnt M ~ c p/q 



We ca n m ake a link with the Brjuno function. The remark in footnote [l] with 



Theorem 1.1 imply that there exists a neighbourhood of the origin where B (lo) 
liie C rit (^) is uniformly bounded. If we are able to prove the previous conjecture 
the same bound will hold in the neighbourhood of every rational, that is we could 
prove B (lo) + In e cr j t (lo) is uniformly bounded for all frequencies. 

In this paper we were interested in finding the optimal scaling law, namely the 



optimality of the exponent 1/q in (4.1), neglecting the behaviour of the proportion- 
ality constants c p / q as function of p/q. It will be interesting to investigate (at least 
numerically) this function, in particular to test its continuity, and to search if the 
function e cr it (Alo) / Alo? is continuous for all rational numbers. 

Appendix A. Numerical precision of the results. 

Here we analyse more precisely the numerical results, pointing out some impor- 
tant remarks concerning the precision of the calculations. 

The destruction of invariant curve is identified by the nonregularity of the fre- 
quency map. In general, this nonregularity is clearly visible (fig 3), and the destruc- 
tion of invariant curves is determined with no ambiguity. But when one searches for 
very fine details, which correspond to small values of the perturbing coefficient e, 
one needs to be able to distinguish between nonregularity of the frequency map and 
the apparent nonregularity resulting from numerical errors in the computations. 

One way to work around this problem is to search always for a regular part of the 
frequency curve. The existence of this regular part is the indication that we have 
not yet reached the threshold for numerical roundoff error. An illustration is given 
in figure ^, where the frequency curve is plotted for e = 0.1, T = 90000, for both 
double and quadruple precision computations. It is clear that with double precision 
computations, we are no longer able to distinguish precisely between regular and 
non regular behavior of the frequency curve, while it is still possible with quadruple 
precision artihmetics. 

As we know that on the regular curves, the FMA algorithm converges towards 



the true frequency as 1/T 4 with the Hanning window of order 1 (p = 1 in Eq.3.5), 
we can check the errors of the FMA method by increasing T. If the variation of 
frequencies increases instead of decreasing, and still present a regular behavior as in 
Figure H.a, it should come from numerical noise, and needs to be checked further. 



A.l. False interpretation of a first result. In fact, the numerical accuracy for 
the iteration of the map is very important when one searches for very small values 
of e, and not taking it properly into account can lead to false results. Indeed, in 
figures 0, we show sections x — xhfp of the Standard Map with e = 0.1 for double 
precision (Fig. fjj.a) and quadruple precision (Fig. |7|.b) calculation. One can note 
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Figure 6. Frequency curve calculated in double precision (a) and 
quadruple precision (b). Section x = 4.408 for the Standard Map 
for e = 0.1, integration time T = 90000. The regular set I reg is 
represented. 



a very large change in the y-range, which implies that the values of y p j q and Vp/ q 
obtained using double precision are largely over-estimated. In this case, we must 
use at least the quadruple calculation, while for the final results of the present 
paper, we had to use the 100-digits multi-precision using the Fortran 90 package 
mpfun. 
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Figure 7. Double precision calculations (a) and Quadruple pre- 
cision calculations (b) for the section x = xhfp for the Standard 
map with e = 0.1. 



Because of the finite precision arithmetic there is a threshold Xthr such that for 
all x below the threshold then sinx = x. When we start with initial condition 
(xoiUo) = o(xthr) then instead of iterating the standard map (1.1) we iterate the 
linear map T l t tn 
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The iteration of T e produces a strange result which is misleading. Using data 
below the threshold value Xthr, the iteration of T e ' m and then the numerical analysis 
of the signal gives the following scaling law: There exists < e* < 1 such that: 

1. if e cr it (uj) > e* then: \ne cr u (Aujr ) ~ - lnAwz + c\, 

\ <i J ( i 'i 

2. if € cr i t (uj) < e* then: lne cr , t I Alue ) ~ - In Awe + C2\ 

\ q J Q q 

namely it seems to be a discontinuity in the scaling law as Figure M shows. 




-7.5 -7 -6.5 -6 -55 -5 -4.5 -4- -3.5 

ln(Aw 0/1 ) 



Figure 8. Change of slope phenomenon for the resonance 0/1. 
lne cr i t is plotted versus lnAwg/i using quadruple precision. 
The linear fit lne C rit = 1.948 In Auq/i + 6.863 and ]ne cr u = 
0.92 In Awo/i + 0.998 are also plotted. 
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Looking at JjJ and comparing with the results of |3(| at first time, using the double 
precision, we interpreted the result of Figure || as follows: 

The scaling law at resonances for the Standard Map is the one pre- 
sented in [Q . The different result of |j36f can be justified saying that 
they didn't go deep enough inside the resonant region and the non 
linearity of the model added some distortion . We are in presence 
of two behaviours: close to the resonant region we find the proper 
law, far away the distorted one. 

In studying further the reasons for this transition, we convinced ourselves that 
this interpretation was not correct. We understood that this behaviour was due 
to the finite precision arithmetic of the iteration of the map. We also observed 
that for values of e smaller than a critical value e^, there was a saturation of the 
measured value of Ay /i = y + (e) — y~ (e), which was also the sign of nu mer ical 
problems (Fig. ^). This was confirmed by the analytic result of Theorem 1.1 for 
the resonance 0/1. 



A. 2. New interpretation. With Lemma 2.3 we proved that the rotation num- 
ber of the standard map is very well approximated by the pendulum frequency. 
Rewriting (|2.6|) as follows 



(A.l) lne cHi ~ 21nAcj 0/1 +C + 21nln 8 /^ 

Ayo/i 

and using the previous remark on Aj/q/Ij it follows that the term lnln : 



12/0/1 



(A.l) is almost constant in the interval of variation of e, say for example e € 
[0.001, 0.01]. A similar result holds for every resonance p/q; then for a fixed numer- 
ical precision we obtain the following scaling law 

2 

lne cri t ~ - In Aw„/„ + c 
q 

To prove that this change of slope is "not real" we increased the precision of the 
calculations (i.e. passing from double precision to quadruple and then to multi- 
precision-100), obtaining that the threshold e* decreases (see Figure 10), showing 
that in the limit of exact arithmetic this threshold doesn't exist. 



Figure pr| shows At/o/i as a function of e C rit for different machine precision. Also 
in this case it is clear that the threshold value under which y is nearly constant 
goes to zero when we increase the numerical precision. 



For the quadruple precision the threshold is Xthr ~ 10 -12 and it goes down to 
order 10 -32 using the multi-precision package 100-digits. This value is sufficiently 
small to do precise calculations in reasonable CPU-times. In fact the threshold 
value for the multi-precision package 1000-digits is order 10 -333 but the CPU-time 
for the calculations increases several order of magnitude, and this was not necessary 
to get good confidence with our present results. 
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Figure 9. In Ay /i as function of e cr it for the fundamental reso- 
nance using the quadruple precision. It is evident that for e < 0.01 
then lnAyo/! ~ —35. 
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Figure 11. Plot of e crit vs lnAt/ /i for different precisions. 
The solid curve is the fitting of the multiprccision-100 data: 
lnAwn/i = a+ -J=, with a = 6.177 and b = -4.804. To 



